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Main Theme We show that in the collapse of a 
Bose-Einstein condensate (BEC) 1 certain processes in- 
volved and mechanisms at work share a common origin 
with corresponding quantum field processes in the early 
universe such as particle creation, structure formation 
and spinodal instability. Phenomena associated with 
the controlled BEC collapse observed in the experi- 
ment of Donley et al [2] (they call it 'Bose-Nova', see 
also 0) such as the appearance of bursts and jets can 
be explained as a consequence of the squeezing and 
amplification of quantum fluctuations above the con- 
densate by the dynamics of the condensate. Using the 
physical insight gained in depicting these cosmological 
processes, our analysis of the changing amplitude and 
particle contents of quantum excitations in these BEC 
dynamics provides excellent quantitative fits with the 
experimental data on the scaling behavior of the col- 
lapse time and the amount of particles emitted in the 
jets. Because of the coherence properties of BEC and 
the high degree of control and measurement precision 
in atomic and optical systems, we see great potential 
in the design of tabletop experiments for testing out 
general ideas and specific (quantum field) processes in 
the early universe, thus opening up the possibility for 
implementing 'laboratory cosmology'. 2 

Theoretical Cosmology in relation to Gen- 
eral Relativity, Particle and Condensed Matter 
Physics For the last half a century the study of theo- 
retical cosmology and high energy astrophysics has re- 
lied largely on general relativity and particle physics, 
while modern and contemporary cosmological exper- 
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2 This essay has the same content as v2 of with a few 
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iments are the fuse and the fuel of these activities. 
We witness the inception of nuclear astrophysics in 
the 50 's, leading to the highly successful theories of 
neutron stars, and particle astrophysics in the 60's ex- 
emplified by the highly successful theory of nucleosyn- 
thesis which helped to establish the standard model 
in cosmology. Establishment of quantum field theory 
in curved spacetime j2] in the 70's laid the founda- 
tion for the study of quantum field processes in strong 
gravity, such as cosmological particle creation in the 
very early universe [7] and Hawking radiation in black 
holes JU (Note these two processes contain very differ- 
ent physics). This pushed the frontiers of theoretical 
inquiries in leaps, back to the period after the Planck 
time. The inflationary cosmology |Hj of the 80's also 
ushered in ideas of particle physics and quantum field 
theory, such as the decay of the false vacuum and the 
capability of vacuum energies driving the universe into 
ultrafast expansion, with scenarios radically different 
from the standard model. 

This side of the story on the progress of modern cos- 
mology with the help of gravitation theory and particle 
physics is well-known. What is perhaps lesser known or 
appreciated is the importance of ideas and techniques 
from condense matter physics and statistical mechan- 
ics in the study of cosmology of the early universe. 
We have seen the relevance of statistical mechanics, ki- 
netic theory, stochastic processes and many-body dy- 
namics in classical astronomy and physical cosmology 
(see, e.g., |1()|). Here we want to emphasize the impor- 
tance of ideas from condensed matter physics in con- 
junction with quantum field theory for treating early 
universe quantum processes, which is believed to have 
played a fundamental role in determining how space- 
time and matter existing in different forms and states 
interplay, transform and evolve. The importance of 
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viewing cosmology in the light of condensed matter 
physics, in terms of taking the correct viewpoints to 
ask the right questions, and approaches to understand 
the processes, has been called to our attention a long 
time ago (see, e.g., ^2])- There were also proposals to 
study cosmological defect formation in helium experi- 
ments and to view cosmology as a critical phenomenon 
|12l Similar efforts aim to identify analogs of full 
cosmological models ^2]. A recent monograph is de- 
voted to the unity of forces at work in He 3 droplets 
jlftj . It should also be mentioned that the proposal of 
sonic black holes II Tj was perhaps the first analog 
model in black hole physics which stimulated recent 
activities in finding similar processes in fluids and con- 
densed matter systems in the so-called analog gravity 
program |T5] . 

Laboratory Cosmology Here we propose using the 
Bose-Einstein Condensate (BEC) and its dynamics as 
another useful venue to 'observe' and probe into some 
fundamental cosmological processes in the early uni- 
verse. Specifically, we analyze the experiment per- 
formed by Donley et al. on the controlled col- 
lapse of a BEC and identify the processes and mecha- 
nisms at work which are responsible for vacuum par- 
ticle creation structure formation and spin- 
odal instability (quenching) in phase transition [23j in 
the early universe. The collapsing BEC is the time- 
reverse scenario of an expanding universe and the con- 
densate plays a similar role as the vacuum in quan- 
tum field theory in curved spacetime. One can under- 
stand the production of atoms in the form of jets and 
bursts as the result of parametric amplification of vac- 
uum fluctuations by the condensate dynamics. This is 
the same mechanism as cosmological particle creation 
from the vacuum, which is believed to be copious near 
the Planck time and during preheating after inflation 
|21j . Some basic ideas common to cosmological theo- 
ries like "modes freeze when they grow outside of the 
horizon" can be used to explain the special behavior 
of jets and bursts ejected from the collapsing BEC. Fi- 
nally the waiting time before a BEC starts to collapse 
obeys a scaling rule which can be derived from simple 
principles of spinodal instability in critical phenomena. 
These examples clearly indicate the great potential of 
a new field of research which we may call "laboratory 
cosmology" , with tabletop experiments designed to test 
the workings of specific physical mechanisms in specific 
cosmological processes. 



BEC Collapse Experiments In the experiment de- 
scribed by Donley et al. 0, a Bose-Einstein conden- 
sate (BEC) in a cold (3nK) gas of Rubidium atoms is 
rendered unstable by a sudden inversion of the sign of 
the interaction between atoms. This is done by alter- 
ing the binding energy at Feshbach resonance with an 
external magnetic field. After a waiting time t co u apse , 
the condensate implodes, and a fraction of the conden- 
sate atoms are seen to oscillate within the magnetic 
trap which contains the gas. These atoms are said to 
belong to a 'burst'. After a time r evo i ve the interaction 
is suddenly turned off. For a certain range of values 
of T evo i ve , new emissions of atoms from the condensate 
are observed. They are called 'jets'. Jets are distinct 
from bursts: they are colder, weaker, and have a char- 
acteristic disk-like shape. 3 

The Model The model is based on the Hamiltonian 
operator for N interacting atoms with mass M in a 
trap potential V (r) = {lu 2 z 2 + cu 2 p 2 )/2, with radial 
p and longitudinal z coordinates measured in units 
4 of a^o, where a^ is a characteristic length of the 
trap, with associated (dimensionless) frequencies lo z — 
^ 'axial i 'w ~ 1/2 and uo p = Uradial/u ~ V2. The inter- 
action is assumed to be short ranged. We introduce a 

3 We call attention to the distinction between the 'Bosc-Nova' 
|21 experiment studied here and other BEC collapse experiments 
[221 1231 . At magnetic fields around 160G, where the effective 
scattering length is of the order of 500ao (and positive) (ao = 
0.529 10 -10 m is the Bohr radius) it is possible to observe 
oscillations between the usual atomic condensate and the molec- 
ular state 1231 with a frequency of oscillations of hundreds of 
KHz 1241 . By contrast, in the 'Bose-Nova' experiment 2 typical 
fields were around 167G, the scattering length was only tens of 
Bohr radii (and negative) and the frequency of atom - molecule 
oscillations may be estimated as well over ten MHz 1251 . While 
coherent resonance between the atoms and the molecules is ex- 
pected to exist for all of these experiments, and has been shown 
to play an important role in the outcomes of some 1251 . we deem 
it unlikely that it plays a dominant role in this experiment other 
than renormalizing the scattering length (For details, see 151). 
Indeed no oscillations are reported in the original experimen- 
tal paper. Instead, as this note shows, the primary mechanism 
for the Bose Nova phenomena is the parametric amplification of 
quantum fluctuations by the condensate dynamics, resulting in 
bursts and jets as particle production from (the squeezing of) the 
vacuum. Recent numerical simulations 1261 and rigorous theo- 
retical investigations 1271 indicating the inadequacy of mean field 
theory seem to corroborate this view. 

4 We use a sign convention such that the effective coupling 
constant is positive for an attractive interaction, and a sys- 
tem of units where the length a; lo , time ti lo and energy scale 
Eho = ^ = M^Q-ho are defined with reference to the average 
frequency ut. We work with units such that these three scales 
take the value 1. 



2 



dimensionless field operator * (r) = a h ^ 2 ^> {x), and 

a dimensionless coupling constant u — \ Tiuja\A U = 
Air (a/a ho ). 

obeys the equation of motion = i H,^ 

and satisfies the equal time commutation relations 
[* (t, r) , 4^ (t, r')] = <S (3) (r - r') . We decompose the 
Heisenberg operator ^> = Q(r,t) + ip(r,t) into a c- 
number condensate amplitude $ and a q-number non- 
condensate amplitude ip, consisting of the fluctuations 
or excitations. 

One crucial point of our analysis is that we shall 
focus on the evolution of the fluctuations for a given 
evolution of the condensate (as extracted from the ex- 
periments); in field theory terms, we shall work within 
the test field approximation. It is fair to say that a full 
theoretical account of the Bose - Nova experiment, de- 
scribing the evolution of both condensate and fluctua- 
tions, does not exist. All we can say with any certainty 
is that we face here a strong back reaction regime, be- 
yond the Hartree - Fock - Bogoliubov approximation 

M 

In contradistinction, the dynamics of the fluctua- 
tions alone may be described, at least in the early 
stages of the experiment, by a simple Bogoliubov ap- 
proximation |29| . We obtain the equation of motion 
for the fluctuation field by subtracting from the full 
Heisenberg equation the Gross - Pitaievsky equation 
(GPE) for <f>. We next parametrize the wave functions 
as $ = <E>oe -ie , ip = i/'oe -40 , where $o and are 
real. During the early stages of evolution, we may re- 
gard the condensate density as time independent, and 
the condensate phase as homogeneous, $0 = $0 i r ) , 
= (t) . We may then write the equation for the 
fluctuation field 



t^--H + E 
at 



o(Vo + ^)=0 (1) 



where Eq = k (oj z + 2to p ), As it is well known, this 
approximation is both "gapless" and "conserving" [301 ■ 
To solve equation we decompose -00 into a self- 
adjoint and an anti-adjoint part ipo = £ + irj, each part 
satisfying an equation 
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5 The squeezing of quantum unstable modes and its back re- 
actions on the condensate has been considered before, e.g., as 
a damping mechanism for coherent condensate oscillations 1311 . 
and applied to the collapse of a homogeneous condensate in |32l . 



^ + [H - E - 2u<P 2 ] £ = 0. 



(3) 



Since the trap Hamiltonian is time - independent, we 
have 



dt 2 



[H-E ]H eff Z = 0. 



(4) 



Here H, 



H — Eq — 2u"I>q. To have an unstable con- 



densate it is necessary that at least one of the eigen- 
values of H e ff is negative; the boundary of stability 
occurs when the lowest eigenvalue is exactly zero. 

One further consideration is that we are interested 
in the part of the fluctuation field which remains or- 
thogonal to the condensate. In the full theory the con- 
densate is the eigenfunction of the one - body density 
matrix with the largest (macroscopic) eigenvalue, and 
the non-condensate is built out of the other eigenstates 
|33| . Since the one-body density matrix is Hermitian, 
they must be orthogonal. The ground state of H e ff is 
certainly not orthogonal to the condensate, since nei- 
ther have nodes. Observe that within our approxima- 
tions, the noncondensate wave function is equivalent, 
up to a normalization, to the phonon operator in the 
particle- conserving formalism |34j . 

If we adopt the values ui z = 1/2, lu p = y/2, relevant 
to the JILA experiment, then instability occurs when 
k = Noa cr it/ah = 0.51. This result compares remark- 
ably well with the experimental value k = 0.55 [2"ll25|. 
as well as with the theoretical estimate presented in 
Ref. |35j . This agreement may be seen as natural, as 
the equations we postulate for the fluctuations may be 
obtained from the linearization of the GPE. In both 
calculations, the geometry of the trap plays a funda- 
mental role. 



Scaling of t co u apse and Critical Dynamics As we 

have already noted, even for condensate densities above 
the stability limit, no particles are lost from the con- 
densate during a waiting time t co a apse . Experimentally, 
tcoiiapse is seen to get very large when the threshold 
of stability is approached from above, in a way which 
closely resembles the critical slowing down near the 
transition point characteristic of critical dynamics. In 
our problem, the quantity which plays the role of relax- 
ation time is the characteristic time e _1 of exponential 
growth for the first unstable mode. This quantity di- 
verges at the stability threshold, which in our analogy 
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corresponds to the critical point. By dimensional anal- 



ysis, we are led to the estimate t co u apse 
to the critical point, we find 



collapse 



tc 



1 



-1/2 



Close 



(5) 



The power law Eq. (|5j) describes with great accuracy 
the way t co u a pse scales with the scattering length; the 
best fit to the experimental data is obtained for t C r%t ~ 
5ms. 
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Figure f: 

In Fig. 1 we plot the scaling law JSJ (full line) derived 
here and compare it with the experimental data for 
Nq = 6000 as reported in Refs. [2] (small black points), 
the tNL ~ {uNq)~ prediction (suitably scaled) as 
given in [32 ESI (dashed line) and the results of nu- 
merical simulations reported in |87| (large grey dots). 
While all three theoretical predictions may be consid- 
ered satisfactory, the tNL ~ (uN ) _1 behavior fails to 
describe the divergence of t co u apse as the critical point 
is approached, and the results of numerical simulations 
reported in j^Zj based on an improved Gross-Pitaevskii 
equation tend to be systematically above the experi- 
mental results, which may be a further indication of 
the quantum origin of this phenomenon |38j. 

We wish to stress that our argument predicts the 
scaling exponent, but not the prefactor; even this ap- 
parently simple aspect of the Bosenova phenomenology 
is surprisingly resilient to theoretical explanation |28j . 
The same scaling law is found from a different perspec- 
tive in 

Bursts and Jets as Amplified Quantum Fluc- 
tuations We now consider the evolution of quantum 
fluctuations, treated as a test field riding on the col- 
lapsing condensate whose dynamics is extracted from 



experiment. The initial state is defined by the condi- 
tion that u = for t < 0; we shall take it to be the 
particle vacuum |0), defined by ipQ (x, 0) |0) = every- 
where. 

One can introduce a mode decomposition of the £ op- 
erator based on the eigenfunctions of [H — Eq] H e f j . 
For short wavelengths A, since H ~ A~ 2 >> 2m<I>q, 
we expect these eigenfunctions will approach the trap 
eigenmodes. The fact that particles in bursts are seen 
to oscillate with the trap frequencies [2] also suggests 
that their dynamics is determined by the trap Hamil- 
tonian. Based on these observations we can assume a 
homogeneous condensate 2u$q ~ k~ 1 o,uj z Nq (t), where 
Nq (t) is the instantaneous total number of particles 
in the condensate. In practice, k _1 is a measure of 
the overlap between the condensate and the excitation 
modes. Therefore, the approximation may be improved 
by adjusting k according to the range of modes where 
it will be applied. 

Let Nq be the initial number of particles in the con- 
densate, and a cr = k/Nq the corresponding critical 
scattering length. Trap eigenfunctions ip^ (r) are la- 
beled by a string of quantum numbers n = (n z , n x ,n y ) . 
The eigenvalues of the trap Hamiltonian are (with 
the zero energy already subtracted) Eft — lu z ti z + 
w p {n x + n y ). There are two kinds of modes, stable 

(oscillatory, or thawed) modes if En > ( uj z , and 

unstable (growing, or frozen) modes if not. In the for- 
mer case we find that, although we assume vacuum 
initial conditions, these modes do not remain empty. 
Up to t collapse, when the number of particles in the 
condensate is constant, the density 



n (r, t) 



(where u>fi 



• 2 ^>i(r)^ (6) 



Ex 



Es 



has a constant 



term and an oscillatory term. This oscillatory term is 
responsible for the appearance of 'bursts' of particles 
oscillating within the trap observed in the Bose-Nova 
experiment 2] . In the WKB limit it describes a swarm 
of particles moving along classical trajectories in the 
trap potential. 

In the opposite case E r i < f^-J w 2, the formulae 
for the density is obtained by the replacement of ua 
in @ by iaa, thus loZ 1 sinu^ - > c^" 1 sinner^. Physi- 
cally their difference is immense. In the first place, the 
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density is growing exponentially, but unlike the previ- 
ous case, there is no oscillatory component, and these 
particles do not oscillate in the trap, in the sense de- 
scribed above. These modes come alive at T evo ive (as 
the scattering length is set to zero), whence they be- 
come ordinary trap modes which oscillate in the trap 
in the same way as the the burst modes . To the ob- 
server, they appear as a new ejection of particles from 
the core of the condensate, which makes up the so- 
called 'jets'. The sudden activation of a frozen mode 
(we are borrowing the language and concept of cosmo- 
logical structure formation) by turning off the particle 
- particle interaction may be described as a "thaw" . 

Observe that in this picture several conspicuous fea- 
tures of jets become obvious. Jets may only appear if 
the turn - off time T evo i ve is earlier than the formation 
time of the remnant. Once the condensate becomes 
stable again, there are no more frozen modes to thaw. 
On the other hand, jets will appear (as observed) for 
Tevoive < ^collapse, when the condensate has not yet 
shed any particles. Also jets must be less energetic 
than bursts, since they are composed of lower modes. 

Beyond t co u apse the number of particles in the con- 
densate, and therefore the instantaneous frequency of 
the excited modes, becomes time dependent. If we con- 
fine ourselves to the early stages of collapse we may 
assume nevertheless that the condensate remains ho- 
mogeneous. Shifting the origin of time to t co u a pse for 
simplicity, we write N (t) = 7V exp (— tjr) (see Fig. 
2). 




Figure 2: 

After expanding in trap eigenmodes we find the two 
kinds of behavior described above. If En > , the 

mode is always oscillatory. If E n < , the mode 

is frozen at t co iiapse, but thaws when exp (— t/r) ~ 
E^a/ auo z . During the frozen period, the modes are am- 
plified, but they only contribute to bursts after thaw- 
ing. If the evolution is interrupted while still frozen, 
they appear as a jet. We therefore conclude that the 
number of particles Nj e t in a jet at time T evo ive is essen- 



tially the total number of particles in all frozen modes 
at that time. This is plotted in Fig 3, for Nq = 16, 000, 
^radial = HO Hz, uo axia i = 42.7 Hz, a = 36a , and 
k = 0.46 , and compared to the corresponding results 
as reported in 
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Figure 3: 

We see that the agreement is excellent at early times 
(up to about 6ms). For later times, this model overes- 
timates the jet population. This is due to the fact that, 
by not considering the shrinking of the condensate, we 
are overestimating the overlap between the condensate 
and the fluctuations, thus delaying the thaw. It nev- 
ertheless reproduces the overall slope of particle num- 
ber with Tevoive t It should also be remembered that 
we are computing the expected number of particles, 
but in the highly squeezed state which results from the 
frozen period, the fluctuations in particle number are 
comparable to the mean number itself. 

It is interesting to observe that it is possible to repro- 
duce the jets within a theory where the GPE equation 
is generalized to include an imaginary three - body re- 
combination loss term We do not regard this as 
necessarily a different explanation, but rather as a dif- 
ferent way of handling the divide between condensate 
and non-condensate (the analogous cosmological prob- 
lem would be whether to consider a spin two fluctuation 
generated during Inflation as a test field on the cosmo- 
logical background, or as part of the geometry). Still 
agreement with the observed jets is obtained only for 
certain ranges of parameters, and it is unclear whether 
any single parameter set gives a satisfactory simulta- 
neous account of all aspects of the experiment |28| . 

In this talk, we have presented a new viewpoint to- 
wards understanding the salient features in the physics 
of controlled collapse of a Bose-Einstein condensate de- 
scribed in the experiment of j2j, i.e, in terms of quan- 
tum vacuum fluctuations parametrically amplified by 
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the condensate dynamics. Even under a number of 
simplifying assumptions, our approach yields results in 
excellent agreement with experiment, particularly in 
the scaling of the waiting time t co u apse and the num- 
ber of particles in a jet. A background field separation 
is assumed here (even though the measured conden- 
sate dynamics contains the backreaction of nonconden- 
sates) so that one can treat these as test-field processes. 
A theoretical treatment of the fully self-consistent dy- 
namics of both the condensate and its quantum fluctu- 
ations during collapse beyond the conventional Hartree 
- Fock - Bogolubov theory remains a worthy challenge 
(see, e.g., an d references therein). 

If our explanation of the salient features of this ex- 
periment is correct one can think of using this process 
to create coherent atoms in highly squeezed states 021 ■ 
This is because the underlying mechanism of paramet- 
ric amplification produces particles from vacuum fluc- 
tuations in squeezed states |43| . 

Our way of thinking here is influenced by insights 
from the quantum field theory of particle creation and 
structure formation in cosmological spacetimes as well 
as theories of spinodal instability in phase transitions. 
One can conceivably design experiments with BEC dy- 
namics to test out certain basic mechanisms and spe- 
cific features of quantum processes in the early uni- 
verse, thus opening a new venue for performing 'labo- 
ratory cosmology'. 
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